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Using a combination of the replica-exchange Monte Carlo algorithm and the multicanonical 
method, we investigate the influence of bending stiffness on the conformational phases of a bead-stick 
homopolymer model and present the pseudo-phase diagram for the complete range of semi-flexible 
polymers, from flexible to stiff. Although a simple model, we observe a rich variety of conformational 
phases, reminiscent of conformations observed for synthetic polymers or biopolymers. Changing the 
internal bending stiffness, the model exhibits different pseudo phases like bent, hairpin or toroidal. 
In particular, we find thermodynamically stable knots and transitions into these knotted phases with 
a clear phase coexistence, but almost no change in the mean total energy. 
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Since the first simulation of knotted polymers in 
1975 [1], the occurrence and behavior of knots in polymers 
has been studied in various contexts. Scanning through 
protein data bases has revealed that several proteins form 
knots [2-4]. In particular, Virnau et al. [5] have reviewed 
the whole Protein Data Bank (http://www.pdb.org [6]) 
and identified 36 different proteins forming relatively sim¬ 
ple knots, none of which features more than five crossings - 
somehow evolution tries to avoid knotted proteins [7]. On 
the other hand, flexible polymers form much more compli¬ 
cated knots, which occur by chance in the swollen [8-10] 
and in the collapsed phase [10, 11]. Essential for the 
existence of these two phases is the excluded volume and 
attraction of the monomers. Already lattice polymer 
simulations [12, 13] show that models integrating self¬ 
avoidance and attraction exhibit a swollen, a globular, 
and a frozen phase. In this work we go a step ahead 
and investigate the knottedness of semi-flexible bead-stick 
polymers. There exists already some works concerning the 
more complex phase space of polymer models that incorpo¬ 
rate bending stiffness [14]. The most comprehensive study 
uses a model for semi-flexible polymers of bead-spring 
type with finitely extensible nonlinear elastic (FENE) 
bonds [15]. By varying the bending stiffness, these mod¬ 
els are able to mimic a large class of polymers, exhibiting, 
for instance, bent, hairpin or toroidal conformations. 

Nevertheless, none of the former work considered the 
knottedness of the polymer over the full bending stiffness 
range, which we will discuss in this Letter. By measuring 
the knot type we found pseudo phases with thermody¬ 
namically stable knotted polymers. The knot type will be 
shown to be an ideal topological order parameter to iden¬ 
tify the knotting transition and, moreover, the behavior 
at the transition from an unknotted to a knotted pseudo 
phase turns out to be surprisingly different from all other 
pseudo-phase transitions of the bead-stick polymer in that 
it does not fit into the common classification scheme of 
first- and second-order phase transitions. 

To model a coarse-grained polymer with an adjustable 
stiffness, we use a modified version of the bead-stick model 
of Refs. [16-18], which neglects the hydrophobicity of 


the monomers. The model then consists of N identical 
monomers connected by bonds with length one. Non- 
adjacent monomers interact via the Lennard-Jones poten¬ 
tial 



where is the distance between two monomers. The 
parameters e and cr are set to one for the rest of this work, 
i.e., energies are measured in units of e and lengths in units 
of cr. The stiffness is modeled through the cosine potential 
adopted from the well-known discretized worm-like chain 
model [19] and defined by 

N-2 

Ebend = E “ COs6>j) , (2) 

i=l 

where Oi represents the angle between adjacent bonds. 
The complete Hamiltonian is then given hy E = Epj + 
/^Ebend? where n parametrizes the strength of the bending 
term compared to the Lennard-Jones potential. 

Topological barriers between the knotted phases forced 
us to apply relatively complex Monte Carlo algorithms in 
order to obtain reliable results. To simulate the system 
in the complete (T, /^)-plane, we used two complemen¬ 
tary Monte Carlo algorithms. The first is a combina¬ 
tion of the parallel multicanonical method [20, 21] and 
a one-dimensional replica exchange [22] in the k, direc¬ 
tion (PMUCA+RE). The second is a two-dimensional 
version of the replica-exchange method (2D-RE), which 
simulates the system in T and k direction in parallel. By 
means of a two-dimensional weighted histogram analy¬ 
sis algorithm [23], we are able to calculate the canonical 
mean values for every point in the (T, /^)-plane within 
the simulated parameter ranges. To generate well equi¬ 
librated results it was necessary to apply quite intricate 
bridge-end and double-bridging moves that respect the 
fixed bond-length constraint, besides the common crank¬ 
shaft, spherical-rotation and pivot moves. The results 
of both methods, PMUCA+RE and 2D-RE, are in good 
agreement with each other. 
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To determine the structural phases we measured the 
total energy (E), both sub-energies (^lj) and (^bend), 
the squared end-to-end distance (i^ee)^ squared ra¬ 
dius of gyration and the eigenvalues of the gyration 

tensor (Ai), (A 2 ) , (A 3 ). The derivatives with respect to 
temperature of these observables mark the locations of 
the different pseudo-phase transitions, which should not 
be mistaken with phase transitions in the thermodynamic 
limit {N 00 ). The different observables give slightly 
different results for finite systems, so that the transitions 
are smeared out. Additionally, we performed microcanon- 
ical analyses as a complementary approach to identify 
the different pseudo-phase transitions. These are based 
on the microcanonical entropy, S{E) = kB lnf2(^), with 
fl{E) being the density of states. In the microcanonical 
framework a peak in the derivative of the microcanon¬ 
ical inverse temperature dPrnicro{E)/dE^ which is itself 
the derivative of the entropy, /dmicro(^) = dS{E)/dE^ 
corresponds to a pseudo-phase transition. A detailed 
description of the framework can be found in Ref. [24]. 
Although it is expected that all but the collapse and 
freezing transitions vanish in the thermodynamic limit, 
they still determine the structural behavior of short and 
midrange polymers and are important for understanding 
the behavior of mesoscopic systems. 

The stiffness-dependent behavior of the bead-stick 
model [Eqs. (1), (2)] for chains with = 14 and 28 
monomers is summed up in the pseudo-phase diagrams 
shown in Fig. 1, which are constructed from the surface 
plots of all measured thermal derivatives and the results 
of the microcanonical analysis. The black lines in Fig. 1 
mark the thermally most active regions and represent 
the location of the pseudo-phase transitions. For high 
temperatures, where the system is entropy dominated 
and the polymer resembles a discretized worm-like chain, 
the conformations are either gaseous-like and extended 
(F) or rod-like (R). When the temperature is lowered, the 
Lennard-Jones energy becomes more important and the 
polymer collapses. One can clearly distinguish two differ¬ 
ent regimes. For small n the polymer behaves similarly to 
a flexible chain with collapsed (C) and frozen (F) pseudo 
phases. This contrasts with larger b:, where the poly¬ 
mer undergoes a first-order-like transition from the un¬ 
structured state (R) to states with differently structured 
motifs (DA/", H, KC^)- For clarity, we omitted in Fig. 1 
some sub-phases, where the shape parameters and the 
microcanonical analysis suggest additional pseudo-phase 
transitions between differently shaped conformations of 
the same motif. Fspecially, the frozen phase F subsumes 
many different crystal-like phases which differ only in 
minor aspects. 

To identify the knotted phases, we measured the knot 
type of the polymer, which turns out to be an ideal order 
parameter. In principle, the knot type, denoted by (7^, 
defines which smooth closed curves can be transformed 
into each other by applying multiple Reidemeister moves. 
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(a) 



FIG. 1. (color online) The black lines sum up all signals 
of the thermal derivative of all measured observables for (a) 
14 monomers and (b) 28 monomers. The pseudo phases are 
labeled as follows: E - elongated, R - rod-like, C - collapsed 
(AC is an artifact of the small chain length, both termini of 
the polymer are aligned), F - frozen, KCn - knotted phase 
with the corresponding knot type, T>N - [N — 1) times bent 
polymers, H - hairpin. We omitted some sub-phases of DA/", 
K 8 i 9 , and F, which do not change the overall picture. 


Practically, this means two knots are not of the same 
type if one of them cannot be deformed into the other 
without cutting and rejoining the curve. The integer 
number C counts the minimal number of crossings for any 
projection on a plane and the subscript n distinguishes 
topologically different knots with the same number of 
crossings. A detailed exposition of mathematical knot 
theory can be found in Ref. [25]. Of course, an open 
polymer cannot satisfy the mathematical definition of 
a knot, unless the termini are closed virtually. In this 
work we applied a closure which is inspired by tying a real 
knot. Both termini are connected by a straight line, which 
is enlarged in both directions to a point far outside the 
polymer. These new termini are then connected through a 
point which lies on a line perpendicular to the connecting 
line of the termini and which is also located far away 
from the center of mass of the polymer. We have tested 
our results with different other closures and obtained 
almost identical results. To calculate Cn, we employed 
a technique described in Ref. [26], which is based on a 
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FIG. 2. (color online) (a) Surface plot of Ap(—1.1) over 
the complete (T, /c)-plane for a 14mer, blue circles corre¬ 
spond to maxima and red squares to minima of ^Ap(—1.1). 
The blue regime, Ap(—1.1) = 9.0546, marks the K3i knot 
phase of the polymer, (b) Temperature profile of {E/N), 
(RI), and Ap(—1.1) at = 1.00 (blue solid lines) and k. = 
3.00 (red dashed lines). In contrast to {E/N) and (RI), 
Ap(—1.1) differs significantly between the knotted {k, = 3.00) 
and unknotted {k = 1.00) phase. 

variant of the Alexander polynomial A(t), 

Aj,{t) = \A{t)xA{l/t)\, (3) 

evaluated at t = —1.1. Ap{t) inherits the ability to distin¬ 
guish different knot types from the Alexander polynomial. 
Strictly speaking, the Alexander polynomial and likewise 
Ap{t) are not unique invariants. However, for simpler 
knots both generate unique results, and they are able to 
distinguish all knots found in this work. 

In the knotted pseudo phases, marked with K in Fig. 1, 
the probability to find the corresponding knot is almost 
one. This means every measured polymer is of that knot 
type and implies that the knots are thermodynamically 
stable. We want to emphasize that this is different from 
the formerly investigated knots which form by chance in 
the swollen or collapsed phase. This makes Ap{—l.l) a 
perfect topological order parameter to distinguish knotted 
and unknotted phases, as is demonstrated in Fig. 2(a). 
In Fig. 2(b) one can see that the behavior of the mean 
total energy (E) and the squared radius of gyration {R‘^) 


is qualitatively very different at the knotting transition 
AC K3i {k, = 3.00) and the freezing transition C F2 
{k = 1.00). On the other hand, the knot parameter 
Ap{—l.l) clearly signals the pseudo-phase transition and 
goes from 1 (unknotted polymer) to 9.0546 (3i knot) only 
in the case of a transition into a stable knot. As expected, 
the 28mer (and even more the 42mer [27]) exhibits a 
richer phase diagram with more complicated knot types, 
see Fig. 1(b). However, the qualitative behavior at the 
phase boundaries turned out to be very similar so that 
we will focus in the following on the 14mer. 

The knotting transitions from one structured state to 
another (e.g., K3i D3) are quite interesting. At first 
glance, one could assume that they behave first-order-like, 
similar to other solid-solid-like transitions at low temper¬ 
atures. However, the microcanonical analysis indicates 
a second-order-like behavior by a peak in dPrnicro{E)/dE 
which is smaller than zero (not shown here). Likewise, the 
canonical probability distribution p{E) does not exhibit a 
double-peak structure, see Fig. 3(a). On the other hand, 
the two-dimensional energy distribution ^(^’lj, ^bend) 
points to a phase coexistence. In Fig. 3(b) one can clearly 
identify two separate peaks, one corresponding to the 
knotted phase and the other to the unknotted phase. Sur¬ 
prisingly, both phases have almost identical mean total 
energy {E) at the coexistence point, and there is almost no 
signal in the total energy and heat capacity at the transi¬ 
tion, see Fig. 4. We thus observe no latent heat where the 
polymer undergoes the transition into the knotted phase. 
Rather, the Lennard-Jones energy Ej^j and the bending 
energy F^bend are transformed into each other [28]. 

This behavior changes if the polymer enters the knotted 
phase from an unstructured conformation. For example, 
staying on the transition line between K3i D3 for 
= 14 and going to higher temperatures, the two peaks 
in p (F^lj, F^bend) Start to merge until they form a single 
peak at the transition between AC K3i, where the 
phase coexistence vanishes. This behavior is similar for 
all other observed knotting transitions. 

The reason for the missing knotted phases in Ref. [15] 
is not rooted in this intricate behavior with “concealed” 
signals, however, but may rather lie in the choice of the 
bond length and the minimum of the Lennard-Jones 
potential rmin which are identical in [15] but different 
by a factor of 1.22 in this work. We have tested a few 
parametrizations of bead-spring models and always ob¬ 
served thermodynamically stable knots except for the 
special parametrizations used in [15]. The observation 
that the polymer minimizes its total energy in bent confor¬ 
mations by maximizing the number of monomers located 
in the Lennard-Jones minima of other monomers suggests 
the following conjecture: if ^ rmin and the bonds are 
flexible enough, bent conformations are energetically so 
strongly favored that knotted states become unlikely. 

In conclusion, we have investigated the complete stiff¬ 
ness dependent behavior of a semi-flexible bead-stick ho- 
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(b) 

FIG. 3. (color online) Energy probability distribution of a 
14mer at a point where ^Ap(—1.1) suggests a knotting tran¬ 
sition {k, = 3.9 and T = 0.0483). Plot (a) shows the one- 
dimensional probability distribution p{E). In (b) the energy 
is split into the two sub-energies, which leads to a double peak 
in the two-dimensional probability distribution p(Elj, Ebend)- 
The left peak corresponds to the knotted state K3i and the 
right one to the bent state D3. The one-dimensional probabil¬ 
ity distribution is exactly the projection of p(Elj, Ebend) along 
the line connecting the two peaks, thus the phase coexistence 
is perfectly hidden in p(E). 


mopolymer. Besides the conformations already observed 
in previous works, we found pseudo-phase transitions 
into novel phases which are characterized by thermody¬ 
namically stable knots and showed that Ap(—1.1) is a 
perfect topological order parameter for these transitions. 
Moreover, at these topological transitions we found phase 
coexistence between structured phases but, surprisingly, 
with equal total mean energy in the two phases. The 
missing signal in the heat capacity implies a pseudo-phase 
transition with a clear phase coexistence but without 
latent heat. At the moment, investigating short single 
polymer chains in experiments seems unrealistic, and it 
appears even more unrealistic to observe structural prop¬ 
erties such as knots. However, preparation and detection 
methods for single polymers at surfaces have recently 
made quite an impressive progress [29, 30], so that one 
can now resolve polymers of down to 20 monomers. A 
sensible next step could be hence an investigation of the 
influence of an adsorbing surface on the occurrence of 
stable knotted pseudo phases for a generic semi-flexible 
polymer. Later on, this could be extended to more realis¬ 
tic synthetic polymer models to guide experiments. 
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fully acknowledged. This project has been funded by the 
European Union and the Free State of Saxony through 
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FIG. 4. (color online) Energy {E) and heat capacity Cv of a 
14mer at T = 0.0483 over the complete range of investigated 
K, values. At both knotting transitions, k, = 2.7 : E 2 K3i 
and K — 3.9 : K3i ^ D3, there is nearly no shift in the energy 
and therefore no signal in the heat capacity that is larger than 
the statistical error. 
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